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Multiphoton resonance in a three - level system with nearly degenerate excited states 
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An analytic study is presented of the efficient multiphoton excitation and strong harmonic gener- 
ation in three-level systems specified by a pair of nearly degenerate, strongly dipole-coupled excited 
states. Such systems are physically formed by the three lowest states in, e.g., the hydrogen atom or 
evenly charged homonuclear diatomic molecular ions under reasonably chosen laser intensities. As a 
detailed analytic result, we found that the laser pulse of photon energy 2.05 eV, duration 0.23 ps and 
intensity 5 • 10 -^ is able to produce complete inversion of the initial population in the hydrogen 
atom through the 5-photon excitation. At the same photon energy, the pulse of duration 0.41 ps 



and intensity 3.44 • 10 



was found to produce the same effect in the molecular ion but through 



the 9-photon excitation. We show that the accompanying scattering of light has very rich spectrum 
differing substantially from that of the two-level system. 



PACS numbers: 42.50.Hz, 32.80.Rm, 33.80.Rv, 82.53.Kp 



I. INTRODUCTION 



Typically, a number of optical photons are necessary 
to cover the energy gap between the ground state and the 
first excited state in most atoms and molecules. Multi- 
photon transition is thus the only way of populating the 
excited states in these materials when exposed to laser 
light. In a nonperturbative analysis, Duvall et.al. [l| 
have, however, proved that in a two-level system multi- 
photon excitation is inefficient because of small multipho- 
ton coupling strengths requiring high laser fields produc- 
ing, in turn, large ac Stark shifts and prompting compet- 
itive processes like multiphoton ionization. As shown by 
Gibson [2|, Q , all these disadvantages of the two- level sys- 
tem are cancelled in a specific class of three-level systems, 
where the multiphoton coupling strengths are greatly en- 
hanced, the ac Stark shifts are greatly reduced and the 
ionization is minimized. The specific three-level systems 
of interest include a pair of nearly-degenerate strongly- 
coupled different-parity excited states separated by many 
photons from the ground state. A physical realization of 
such systems are, e.g., all evenly charged homonuclear 
diatomic molecular ions, like iVj for example, as re- 
minded by Gibson. The recent numerical calculations by 
Gibson for a ID model of this ion have shown that its 
ground state is distanced by as much as over 0.66 a.u. 
(over 18 eV) from the excited pair of states at the inter- 
nuclear separation of the order of 3.5 a.u.. At this in- 
ternuclear separation, the splitting of the excited states 
amounts to 0.0167 a.u. (0.454 eV) and the dipole matrix 
element between them is 6 times greater than that be- 
tween the ground state and the lower excited state. By 
numerical calculations Gibson found that, at laser inten- 
sities of the order of 10^^ -^, nearly complete inversion 
of the population could be obtained in the 11-photon and 
12-photon excitation processes. Also, an important find- 



ing was that no other discrete ionic states, besides the 
three lowest ones, were populated during the excitation 
and practically no ionization occurred at intensities not 
exceeding substantially 10"'^^ -^. All this means that, 
up to the mentioned laser intensities, the ion behaves 
like a perfect three-level system. Thanks to the efficient 
multiphoton excitation, this ion was predicted by Gib- 
son to be a source of strong generation of harmonics of 
the incident laser light. At the above intensities, the 
harmonic generation by the ion differs in its mechanism 
from the generation by atoms because in the latter case 
the atomic continuum is strongly engaged through the 
three-step mechanism including ionization, acceleration 
of the ionized electron and, finally, its recombination [4|. 
In the case of atoms, e.g., the hydrogen atom with its 
three lowest states IS', 2P and 25", a substantial decrease 
in the laser intensity by almost two orders of magnitude 
is required to avoid the ionization effects and to approx- 
imate the atom by an effective three-level system. 

In the present paper we give an approximate analytic 
solution to the problem of the above specific three-level 
systems strongly driven by a laser field. Our analytic so- 
lution will be shown to confirm the numerically-based 
predictions of Gibson concerning efficient multiphoton 
excitation and strong harmonic generation. The chain of 
dipole couplings in the system under the present study 
forms a F-type configuration (see FiglT]). This config- 
uration differs from the A-type and ladder-type config- 
urations we have previously studied [5| under the op- 
posite assumption of weak multiphoton excitation and 
by a different method based on a set of two differen- 
tial equations for appropriately defined ratios of the level 
population amplitudes. At weak excitations, this set was 
made decoupled and one equation was transformed into a 
quadratically nonlinear Riccati-type equation. Thus, the 
present paper considerably extends our previous treat- 
ment by including a different system, a different approach 
and higher laser intensities. 
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II. THEORY 



The model 



The system of interest, shown in FigH], includes three 
states \j) [j = 1,2,3) of eigenfrequencies ujj and defined 
parities (either "+" or "-"). The parities of states |1) 
and 1 3) are assumed to be the same and opposite to 
that of state |2). This system, with its initial population 
in |1), is driven by a light of linear polarization along 
the z-axis, carrier frequency wq, electric field amplitude 
Eq and switching on (off) function < J{t) < 1. The 
interaction is taken in the electric-dipole approximation 
V{t) = ~ezE{t), where E{t) = Eofit) coa{ujot). Instead 
of Vij (t) , we introduce the instantaneous Rabi frequencies 

for the 1 — 2 and 2 — 3 couplings, namely V.{t) = ^^ = 

nnfit) cos{uot) and M(i) = -^^M. ^ Mfl/(i) cos(woi), 
where flR = ti^ and Mr = H^jEo .^^.^ ^^^ standard 
time-independent Rabi frequencies while ^jk = {j\ez\k) 
is the dipole matrix element for the j -^ k transition. 
Throughout this paper the duration of the light pulse is 



|2> 



(-/+) 



M(t) 



(+/-) 



(1^3 



-|3> 
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CJ2 



n(t) 



-faji 



Figure 1: The three-level F-type system with nearly degener- 
ate upper states in a laser field of frequency a;32 ^ a;o <^ lj2i- 
il{t) and M{t) are instantaneous Rabi frequencies for the 1-2 
and 2-3 dipole couplings, respectively. The dot points to the 
initial population. 



taken to be much shorter than the excited state lifetimes, 
so the function of the system under interaction is repre- 
sented by ^(i) = J2j=i^ji'^)\^)- From the Schrodinger 
equation for \['(t) one obtains the differential equations 
for the time-dependent population amplitudes bj{t): 



ibi — ujibi — fl{t)b2 , 

ib2 = LU2b2- n{t)bi ~ M{t)b3 , 

ibs = ^363 - M{t)b2 ■ 



(1) 



We perform a series of transformations on this set. 
Four of them are those suggested by Gibson 3] for the 
system with degenerate upper states (0^2 — ^3)- In turn, 
the first transformation is 



bj e*"^* , 



(2) 



the second 








Cl = 


= Cl, 




c± = 


= ^(^^±^3) 


the third 







di = Cl e 



d± 



c± e 



T#(t) 



and the fourth 



di 

^2,3 



dl , 
-i=(d+±d_) 



(3) 



(4) 



(5) 



where cuap = uja—^p, (f'it) — L _„ M{t )dt and <o is the 
turn-on time of the pulse. To 1;hese transformations we 
add one more transformation 



X = di, 




y ^ d2 exp 


JQ 


z = d^ exp 


iu}32 / cos^ {(t){t' ))dt 
Jq 



(6) 



By these transformations the initial set of equations for 
bj is replaced by a new set for x, y and z: 

ix = F{t)y + G{t)z, 
iy = F*{t)x + H{t)z, 



where 

Fit) 

Git) 



iz 

- exp 
'i exp 



G*{t)x + H*it)y., 



(7) 



--(W21 +W3i)i 



'^i^2l+^il)t 



Hit) = -i'^q-\t)sini2m) 



g(t)O(i)cos(0(t)),(8) 

g-i(t)n(t)sin(0(t)),(9) 
(10) 



with 



g(t) = exp 



-^32 / cos(20(i ))dt 



(11) 



Obviously, set d?]) for a;, y and z is as exact as set dD) 
for bj. The advantage of the set for x, y and z is that it 
makes it possible to identify multiphoton resonances in 
the system. These resonances are rooted in the coupling 
terms F(t), G(i) and Hit) being nonlinear in M(i). 

Having found x, y and z we return to the amplitudes 
bj by the inverse transformations 



bi = e-*"i* 



62 
63 






(12) 
e --n-\t)iFit)y + Git)z), (13) 

e— ^*0-i(i)(g-2(t)i^(f)z-hg2(t)G(t)y).(14) 



From the initial conditions for bj{t) (t'i(O) = 1, 62(0) = 
63(0) = 0), one obtains the following initial conditions for 
X, y and z, namely a;(0) = 1, y(0) — z(0) = 0. We will use 
the above amplitudes to calculate the population dynam- 
ics, |fej(i)P, and the coherent part of the correspondence- 
principle spectrum of scattered light. The total spec- 
trum, S'(aj), is defined by [6| 






Ate^-\-^,{t)\ez\^,{t)) 



(15) 



where tp is the duration of the incoming pulse and 
^ j{t) stands for the Schrodinger wave function satisfy- 
ing the condition ^j(O) = \j). The diagonal term in 
Eq. p^ . j = 1, defines the coherent part of the spec- 
trum, Sc{i^)^ and the rest is the incoherent part, Sinc{^)- 
The coherent part results from the average dipole mo- 
ment d{t) = (*i(i)|ez|*i(t)) = 2^ll2Re{bl[t)b2{t)) + 
2 H23 Re {bl{t)b3{t)), while the incoherent part from the 
dipole fiuctuations around the average value. In all our 
calculations we will take many-cycle pulse, i.e., much 
longer than the optical period T = —, meaning the 
pulse shape function f(t) to be slowly varying in time 
when compared with cos(a;oi). This assumption will al- 
low us to remove the slow function f{t) from the inte- 
grand when multiplied by a fast function. Thanks to 
that, fi{t) cos{4>{t)) and il{t) sm{(j){t)) , being the major 
parts of F{t) and G(t), respectively, can be expanded in 
harmonic functions as 



2k 

Jk-i{p) + Jk+i{p) = —Jk{p) ■ (20) 

P 



Exact degeneracy of the excited states (1032 = 0) 



In the case of exact degeneration of the upper levels, 
q{t) — 1 and H{t) — and then set ([7]) reduces to 



a = F{t)y + G{t)z, 
iy = F*{t)x, 
iz — G*{t)x, 



where now 



^R 






n=l,3,5,... 



I i{nujo—ui2l)t I —i{nLOo+uJ2l)t\ 
■" p=2,4,6,... ^ ^ 



' i{pU}Q~U}3l)t _ ^-i(pUo+U3i)t 



(21) 



(22) 



(23) 



n{t) cos{cf){t)) 



ri(t) sin(0(t)) = 



n 



R 



Mn 



2oJo^J2^2k+l)J2k+i 

^ k=Q 

cos [{2k + l)ujot] , 

^-.^f:m.hk(^m 

sin {2kujf)t) . 



cjo 



I{t) 



(16) 



(17) 



The above expansions result from the known relations for 
the Bessel functions Jk{x) J\: 



When Mb. -^ 0, we have F{t) = -f7(i)e-''^2it ^nd 
G{t) = and then equations (|2ip are made the standard 
equations for the two - level system under one-photon 
resonance. 

As seen, F{t) determines the odd - photon excitation of 
|2), while G{t) - the even-photon excitation of |3). Under 
the odd-photon resonance (a;2i = Nljq + Sn, \Sn\ ^ ^O: 
N = 1, 3, 5, • • • ), F{t) splits into the slow part F^ and 
the rapidly oscillating part f^, while G{t) — Gn always 
remains a rapidly varying function. Precisely: 



cos(psina) ~ Jo(/o) + SYJ >^2fe(p) cos(2fca) , (18) 

fc=i 

sin(psina) = 2 ^ J2fc+i(p) sin((2fc + !)«) , (19) 



fc=0 



Mr \ Wo 



-iSj\it 



-ajs^e 



-iSjsft 



(24) 



TN 



Mr 






(25) 



r 



G 



N 



^R _-iS^t 



-"'m' 



p=2,4,6,--- ^ ^ 



' ^i(p-N)uiot _ ^-i{p+N)uiot 



(26) 



In the other case of even - photon resonance (wai = Pujq+ 
Sp, \Sp\ <C Wo, P — 2, 4, 6, • • • ), the functions F{t) and 



G{t) change their roles. Now, G{t) is composed of its slow 
part Gp and the fast part gp, while F(t) — Fp remains 
a rapidly varying function: 



Gp = -P.o^Jp (^m , e 

Mr \ UJo 



-ape 



-iSpt 



(27) 



gp = 



-'wr'- 




pjf 



(Mr 
\ Wo 



m 



^i{p-P)uot 



E 

0=2,4,6, 



PJp ( — fit) 



V ^0 



-^-i{p+P)ujQt 



(28) 



r 



Fp 






E 



nJn f[t) 

\ Wo 



f ^i{n~P)uot I -i{n+P)iJot\ 



(29) 



The slow functions F/v and Gp are seen to have the 
same structure and become constant for exact multipho- 
ton resonance {Sn = Sp = 0) caused by a square pulse 
(/(i) = 1). These constants are the resonant multiphoton 
Rabi frequencies identified by Gibson earlier (Eq. (14) in 
Q). On the other hand, the fast functions oscillate with 
either even multiplicities (/at, gp) or odd multiplicities 
(Gat, Fp) of the carrier frequency wo of the incident light. 
Under a given multiphoton resonance, equations (pi]) are 
transformed into 



(30) 



Hn — 


{Fn + fN)yN + GnZn , 


Wn = 


{Fn + f*N)xN , 


izN = 


G^xn 


photon 


excitation {N = 1, 3, 5, • 


lip = 


Fpyp + (Gp + gp)zp , 


iyp = 


Fpxp , 


izp = 


{G*p+g*p)xp 



(31) 



for the even-photon excitation (P — 2,4,6, •••). Thus 
far, no approximation (besides the assumption of slow 
pulse envelope) has been made, when manipulating the 
equations. 

In the lowest-order (rough) approximation, we are 
tempted to neglect all rapidly oscillating terms in equa- 
tions ([50]) and (|3ip in complete analogy to the standard 
rotating wave approximation (RWA) known from the 
one-photon resonance [sj, |9|. In such an approach, 
equations ((30|) are reduced to Hn — — aAre~*''"*j/Ar, 
Wn = —a Are'*"* a; AT, izff = 0, while equations (|31l) to 
Up = —ape~^^''*zp, izp = —ape^^^^xp, iyp — 0. In 
either case, the first two equations are structurally the 



same as those for the two-level system under RWA, but 
now with the multiphoton couplings uk (K = N, P) 
defined by equations ([M)) and (P7|) . Obviously, in 
the limit Mr — >■ 0, only ai takes a non - zero value 



ttK 



CtBjjt 



2 — ^K.i where Sa.b is the Kronecker symbol) . 



Choosing, for example, a square pulse (/(t) — 1) we find 



XNit) 



XNiO)cosiANt) + ^{^^NiO) 



+ aNyN{0)) sm{ANt) ]e' 



ywW 



yN{0)cos{ANt) 



A 



N 



i^VNiO) 



aNXN{0) ) sin(ylAft) e* 



ZN{t) = ZNiO), 



where An = 



'AT 



(32) 



has the meaning of an ofi'- 
resonant multiphoton (odd - photon) Rabi frequency.lt 
results from the derivation that the solutions (|32|) have 
to change slowly on the scale of the field period, so the re- 
striction of their applicability is ^Af "^ wo. For the even - 
photon resonance, xp is obtained from xn by changing 
N —> P and y -^ z simultaneously; zp is obtained from 
j/A? in the same way. Obviously, yp(t) = yp{0) in the 
latter case. For the initial conditions 6i(0) = 1, 62(0) — 
63(0) = 0, equivalent to x{0) = 1, y{0) = z(0) = 0, we 
combine the solutions ([5^ with the transformation rela- 
tions (|T2|) - HH) and obtain the time-dependent popula- 
tion amplitudes bj . In the case of odd-photon resonance 



(TV = 1,3, 5,. 



^r'w 



bi-\t) 



■ ), these amplitudes are 

cos{ANt) + i-— — sin(^Afi) 

2An 

i— — sin(^Afi) cos(0(i))e*^^ 
An 



J{Nuio + - 



-bJ2)t 



-CJ2)t 



kW 



it) 



a-N 



A 



N 



sm{ANt) sin(0(i))e'(^^-"^)* . (33) 



For the even-photon resonance (P = 2, 4, 6, • • • ), one ob- rapid parts are found to satisfy the set 



^2 from 63 ' and 63 from &2 



tains bi from bl , 
by changing iV — > P in equations (|33|) . At any time t, the 
above amphtudes satisfy the condition of normahzation 
of the total population probability, |&i(i)P + l&2(i)P + 
I ^3 (OP ^ 1- With these amplitudes we can analyze 
not only the population dynamics, |5j(t)p, but the co- 
herent spectrum of scattered light as well, S(u))/uj'^ — 
I /p'' (iie*'^*(i(t)p, where the average dipole moment is now 
reduced to d{t) = (*i(i)|ez|*i(i)) = 2^12^6(6*62) be- 
cause Re (6263) = as results from (|33|) . For A^-photon 
excitation {N = 1,3, 5,---), we obtain from equations 
((33)) and the Fourier-Bessel expansions (TTSl) , (|19l) the fol- 
lowing induced dipole: 



iPx = iJNfNit) + ZNGN{t) . 
i$y = XNfNit), 
iPz = XNG*j^{t) . 



(35) 



Having in mind that Xn, Vn and zn are slow functions 
with respect to Jn and Gat, this set is easily integrated. 
When the result of this integration is substituted to the 
equations for the slow parts and then time-averaging of 
the right-hand sides is made, one obtains 



IXN 

iZN 



(A/„ + Ag„) XAT - QNC 



'iSjyt- 



aNe'-^'^^XN , 



VN , 



(36) 



d(~)(t) 



/^12 



ajv 



2A 



N 



fe=0 ^ ^ 



( -^ [cos(A^ - 2k)ujnt + cos{N + 2k)ujat\ 

L Am 



1- 



JN 



2An 

+ cos{{N + 2k)ujo~'2AN)t 
Sn 



cos {{N - 2k)LJo - 2An) t 



1 



2A 



N 



COS {{N ~ 2k)LJo + 2 An) t 



+ cos((A-|-2fc)a;o + 2Ajv)t] I 



(34) 



where ao — ^ and afc>o = 1, by definition. In the case 
of even-photon excitation (P = 2, 4, 6, • • • ), S^^ (i) is ob- 
tained from the above (i^^^(i) in the following steps: first 
we put 1 in place of all at, then replace N ^ P and 
2fc — >■ 2fc + 1 and, finally, change the sign -I- -^ — at 
the second cosine in each bracket. These average dipoles 
suggest that odd-order harmonics will be revealed in the 
coherent spectrum Sci^^) with each harmonic having a 
triplet structure, in general. At exact multiphoton res- 
onance {5k ~ 0, where K = N,P), the triplets convert 
into doublets with the field-dependent separation Aapc 
between the doublet components. In the limit M^ — ^ 0, 
we recover from Eq. (|M)) the known coherent spectrum, 
i.e. the doublet with separation 2flji around the perfectly 
resonant incoming frequency wq — ll!2i. 

An improved solution to equations (|30|) is found by 
applying the approach of Avetissian et al. [13, lU . Ac- 
cording to this approach, x, y and z are separated into 
their slow and rapid parts, e.g., x{t) — x{t)+j3x{t), where 
x{t) is the time average of x{t) while Px{t) is a rapidly 
oscillating function on the scale of the incident field pe- 
riod T = —. Then, we substitute the two-part expres- 
sions for cc, y and z into (j30p . separate the slow and 
rapid oscillations and neglect small terms Fj^Py ^ P^ 
and F^Px <C /3j, in equations for the rapid parts. The 



where 



= -ifN{t) / fUt'W 



(Sjv|<wo 



2Nuj, 



f^\ 



\MrJ 



E 



n 



n=l,3.5, 



,._^M^m]m 



Ac 



\Sn\^1^0 r,AT. . ( ^R 



^G^^(i) / G%it')dt' 
2 



= 2Nujo 



\Mr 



E 



p=2,4,6, 



p^ ~ N 



Jl f ^fit) ) (38) 



2 p 



Wo 



have the meaning of the dynamic Stark shifts. By the 
substitutions xn — wexp[— i(A/„ -\- A.Qj^)t\ and ijn — 
wexp [iAy„t], the first two equations in (1361) are reduced 
to iii = — a^r exp [— iJjvi] u, iv = — cat exp [i^^vi] w, 
where 5fq — 5n — (2A/„ -I- Ag„) is the Stark-modified 
N - photon detuning. The equations for u, v have the 
same structure as the previously discussed set i xn — 
— Oat exp [— j(5Art] j/AT, iyN — —a^ exp[—iSNt]xN. Thus, 
u and V are obtained from the first two equations (|32l) . 
by replacing xn ^^ u, Vn ~^ v and also Sn -^ Sn, 



A 



N 



A 



N 



'N 



'(N) 



As a result, the popula- 



tion dynamics, 6^- , still results from equations (j33p 

after the replacements Sn ^^ Sn, An — > An- The same 
replacements must be done in equation (1341) for the av- 
erage dipole. With the above Stark - type corrections, 
the applicability range of the solutions is now An <S^ wq • 
Some corrections to both the population dynamics and 
the average dipole are expected to come from the rapidly 
oscillating parts jS^, (3y and /3z oi x, y and z. These rapid 
parts are to be calculated from equations ([55)) after the 



slow parts Xn, Vn and zm have already been found. Ob- 
viously, what has been said above holds for the even - 
photon excitation as well. 



Near degeneracy of the excited states 



is a correction to the coupling parameter oat, while 

|w32|^tJ0 



A 



H 



-iH{t) f H*{t')dt' 



n=l,3,--- 

is an additional Stark shift. 



(40) 



The theory of Section IIIBI can straightforwardly be 
adapted to the case of nearly degenerate upper states, 



^32 



< 1. 



Af- 



i.e., those satisfying the condition 

ter using equation P^ . q{t) of equation pT|) can 
be reduced to the form q{t) — exp(^uJ32t), where 

^32 = i032Jo (^/(i)) = C^32 - 2A(i) with A{t) = 

^(l- Jo (^/(O))- As a resuh, F{t) is stiU given 
by equation (|22p and G{t) by equation ([23|l with the re- 
placements 0-121 -> CJ21 = a;2i + A(t) and W31 — > wai = 
cjsi — A(t), respectively. Obviously, i?(i) of equation pUj) 
differs from zero now but is a rapidly oscillating function 
due to (J19p and our approximation for g(i). In the case 
of multiphoton resonances , Sn ^^ Sn — ^21 — Nujq and 
5p ^Sp ^Lj3i~ PujQ in equations (gl]) - ([25]) and (gT]) - 
([^. respectively. Also, the term H{t)zN {H*{t)yN) 
must be added to the right-hand side of the second (third) 
equation of the set (j30[) . Similar extension must be ap- 
plied to the set (|3T|) but with the change N ^ P now. 

These changes only slightly modify the solutions of Sec- 
tion HlB] For example, making the rough RWA, consist- 
ing in neglecting the rapidly varying terms /jv, Gn and 
iJjv in the extended set pop . one arrives (for f{t) = 1) at 
equations ([32]) . ([33|) and (p4| with the change Sn ^ Sn. 
On the other hand, when applying the improved method 
of Avetissian et al., we again approach the set (j35|) in 
which iPy has the extra term ZMH{t) on its right - hand 
side, while 1/3^ the extra term yNH*{t). Also, we recover 
the set (|36| with Sn —>■ Sn and with some coefficients 
appropriately redefined. Precisely, ojv — >■ o-n + Q^Af in the 
first equation, un — > aw + ctN and A/„ -> A/„ — Ah 
in the second equation, and Ag„ — > Ag„ + A/f in the 
third equation, where 



UN 



iGN{t) I H*{t')dt' 



kj32|<CtJo ^R 



Jp-N 



p=2,4,--- 

/2Mr 



\ i^O 



-fit) 



Mr 

UJQ 
P 



fit) 



p~N 



(39) 



III. APPLICATIONS 

As an example of a three-level system with degenerate 
upper states, we take the three lowest states in the hydro- 
gen atom: |1) = llS*), |2) = |2P) and |3) = \2S). In this 
case, the excitation energy is UJ21 = 0J31 — | a.u., while 
the dipole matrix elements amount to /ii2 = 0.745 a.u. 
and /X23 = —3 a.u., so /J.12/M23 = —0.248. For a square 
pulse, the system is expected to respond according to 
equations ((33)) and (|34)) in which Sn ^ Sn = 5n — 
(2A/„ -I- Agjv), where A/^^ and Ag^^ are to be calcu- 
lated from equations ([57)1 and ([551) . respectively. Thus, 
the N-photon resonance should occur at the light fre- 
quency Wo = ;^('^2i - (2A/„ + Ag„)). We take a field 
allowed by the condition of applicability of equations ([55)) 
and (|34p. For a given N-photon resonance, this condi- 
tion reads \aN/ujo\ = l^^'^N (^) I < 1- Assuming 

N ^ 5, for example, we find the ratio jas/cjol < 0.1 if 
Mr/ojq < 3.5. So, we choose Mr/uq = 1.5 and ob- 
tain, for this field-strength parameter, the following Stark 
shifts A/5 = -0.00926^0 and Aqs = -0.00646 wq, giv- 
ing the 5-photon resonant frequency wq = 0.0754 a.u.. 
The calculated wo leads to the electric field amplitude 
Eq = 1.5a;o/|M23| = 0.0377 a.u., corresponding to the 



light intensity 1.42 • lO"'^ a.u. (4.99 • 10 



ll3 w 



).For the 



above atom-field parameters, we show in Fig. 2 the evolu- 
tions of the level populations and the coherent spectrum 
of the scattered light. Precisely, Fig. 2a and Fig. 2b show 
the evolution of the ground-state population, |6i(i)p, and 
the evolution of the total population in the excited states, 
\b2it)? + |^3(t)P- The bold solid curves in Fig.2a and 
Fig. 2b were obtained from the simple analytic solution 
(|33| . The dotted curve with small fast oscillations in 
Fig. 2 a is the result of applying the improved approach 
of Avetissian et al. when finding the fast parts Pxi (3y 
and Pz of Xn, Vn and zn, respectively (see the remark 
at the end of section Hi Bp . In Fig. 2b, the dotted curve 
with rapid oscillations conies from direct numerical inte- 
gration of equations ([T]) for the field being turned on by 
10 cycles according to sin^(a;oi/40), where a;oi/40 < 7r/2. 
As seen, the approximate analytic evolutions (Fig. 2a) im- 
itate quite well the numerical ones, though some differ- 
ence in the periods of slow oscillations between the an- 
alytic and numerical curves must be noted. An inipor- 
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Figure 2: (Color online) Five-piioton (A^ = 5) IS — )> 2P res- 
onance in the hydrogen atom produced by the laser field of 
strength Eo = 0.0377 a.u. and frequency ujo = 0.0754 a.u. (a) 
and (b), temporal evolutions of level populations, (c) coher- 
ent spectrum of scattered light . Bold solid curves correspond 
to approximate analytic solutions (either Ea.(|33l) or Eq. (|34p V 
The small rapid oscillations superimposed on the bold curve 
in (a) come from the approach of Avetissian et. al. while those 
in (b) come from numerical integration of Eq.([l| for the field 
with 10-cycle turn-on time and then of constant amplitude. 
The turn-on function was sin'^ (cjot/40) with the restriction 
uJot/AQ < ^. In (c), the numerical spectrum was obtained for 
the field with the same turn-on function as in (b). 



tant conclusion resulting from Fig. 2a, b is that the sim- 
ple solution (1551) describes satisfactorily well the main 



Figure 3: (Color online) Nine-photon (A'^ = 9) resonance in 
the molecular ion ylj"*" produced by the laser field of strength 
Eo = 0.099 a.u. and frequency ujo = 0.0754 a.u. {uio the same 
as in Fig[2|. The description as in Fig[2] 



trend observed in the numerical evolution, i.e., the pos- 
sibility of achieving the total inversion of the population 
through niultiphoton (5-photon) resonance. For the pa- 
rameters of Fig. 2, 225-cycle pulse (0.45 ps) turns out to 
be long enough to execute one trip in the population evo- 
lution, i.e., to move all the population from the ground 
state to the excited states and then back to the ground 
state. Under the condition of complete population in- 
version, Fig. 2c shows the coherent part of the spectrum 
of scattered light (solid curve - Ea. (|M)) . dotted curve - 



numerical integration of Eq.([T])). Except the height of 
the background, the analytic and numerical curves agree 
with respect to the number of peaks, their positions and 
relative heights. 

To present the results for a three-level system with 
nearly degenerate excited states we use the system pa- 
rameters calculated by Gibson [3] for his one-dimensional 
model of evenly charged homonuclear diatomic molecu- 
lar ion A2^ , where ^ = N, O, I. It results from figures 
10 and 11 of Gibson [Si] that, at the internuclear separa- 
tion 3.5 a.u., the three lowest levels are so spaced that 
W21 = 0.6685 a.u., 0^32 — 0.0167 a.u. and, moreover, the 
dipole matrix elements are /ii2 = 0.503 a.u. and fi23 — 
3.033 a.u.. Now, Sn ~^5n ^Sn + A~ (2A/„ + Ag„) in 
equations ^ and §^, with A = ^ fl - Jo (^))- 
As an illustration, let us consider the 9-photon resonance 
which is predicted at ojq — |(cj2i + A — (2Ajg -I- Agg)). 
We find that jag/ojol < 0.1 provided that Mr/ujq < 7.15, 
so we choose Mh/ujq — 4. For this value of Mh/ujo, we 
calculate A = 0.0069 a.u., A/g = -0.0155 wo and Agg = 
— 0.014 wq. These values recommend to take the laser fre- 
quency Wo = 0.0754 a.u., i.e., the same as that in the case 
of the 5-photon resonance in the hydrogen atom. The fre- 
quency Wo gives the scaled separation W32/W0 = 0.22 and 
the peak electric field Eq = 0.099 a.u. translating into the 
light intensity 0.0098 a.u. (3.44 • lO" ^). As found by 
Gibson [3j, in fields like that practically no other levels 
of a real molecular ion, besides the three discrete lev- 
els of the model, take part in the interaction. With the 
above ion-light parameters, we show in Fig|3] the popu- 
lation evolutions and light scattering spectra. The quali- 
tative conclusions resulting from Figl3] for the molecular 
ion match those drawn from Fig [5] for the hydrogen atom. 
A quantitative difference between the figures is that, for 
the molecular ion, a longer pulse of 405 cycles (0.82 ps) is 
required to execute one complete trip in the population 
evolution (ground state— >-excited state— aground state). 
The other difference is a richer spectrum of scattered light 
including 12 components instead of 6 components met in 
the case of the hydrogen atom. 

We add Fig H] to show that, by lowering the order of 
resonance to iV = 3, specifically for the hydrogen atom, 
a better agreement was found between the analytically 
and numerically calculated population evolutions. 

We wish to stress that by ignoring the third states in 
the systems under study one obtains quite different re- 
sults. For the hydrogen atom and the same field strength 
and frequency as in Figl^l we found by numerical inte- 
gration that the ground state in the two-level system is 
only weakly depopulated, and the spectrum includes only 
three components (wo, 3 wq, 5wo) with the 3-rd and 5-th 
harmonics being much weaker. Qualitatively the same re- 
sult was found for the AP' ion and the same field strength 
and frequency as in FiglS) The found weak depopulation 
of the ground ionic state resulted in a drastic reduction 
of the number of spectral components in the scattered 
light (to two only) and in lowering the height of the third 
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Figure 4: (Color online) Population evolutions as in Fig[2]but 
for three-photon (A'^ = 3) resonance in the hydrogen atom 
(cjo = 0.1255 a.u., Eo = 0.0314 a.u.). In (b), the numerical 
curve is for field with 2-cycle turn-on time. 



harmonic. In FiglS] we show the numerically calculated 
spectra of scattered light for the above mentioned two- 
level systems obtained from our original three-level sys- 
tem (FiglD) by ignoring the third state of the same parity 
as that of the ground state |1). Thus, the importance of 
the third state is that the three-level system interacts 
much stronger with laser light giving stronger multipho- 
ton excitation and harmonic generation, because of its 
strong coupling to the second state by large dipole mo- 
ment. 



IV. SUMMARY 

In a specific three-level F-type system, we have treated 
analytically both multiphoton excitation and the accom- 
panying generation of harmonics of the incident light. 
The specificity of the system was that it contained a pair 
of perfectly or nearly degenerate excited states separated 
by many photons from the ground state and, moreover, 
the dipole coupling between the excited states greatly ex- 
ceeded that between the ground state and the lower ex- 
cited state. Up to 100% 5-photon excitation in the hydro- 
gen atom and 9-photon excitation in the A2 molecular 
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ion were found by subpicosecond laser pulses of moder- 
ate intensities (5 • 10" -^ (hydrogen) and 3.44 • 10^'^ -^ 
(molecular ion)) at the frequency 0.0754 a.u. (2.05 eV 
photon energy). The calculated coherent spectrum of 
the scattered light was found to be composed of a num- 
ber of odd-order harmonics and it was much richer and 
better pronounced than that from the two-level system 
obtained by rejecting one excited state of the same parity 
as that of the ground state. The analytic results agree, 
qualitatively at least, with those of numerical integration 
of the Schrodinger equation. The quantitative differences 
are that the numerically and analytically calculated pop- 
ulations oscillate slowly with slightly different frequen- 
cies, and the spectra of scattered light though coinciding 
in peaks have different heights of their backgrounds. A 
possible reason for different frequencies of slow oscilla- 
tions is that, in the analytic calculations, we did not find 
sufficiently precisely the frequency cjq generating multi- 
photon resonance. This supposition is maintained by the 
observation that changing slightly the calculated wq we 
were able to achieve better agreement between the ana- 
lytic and numerical curves in figures [5] and |31 Also, by di- 
minishing the multiphoton order of resonance and, thus, 
the required laser intensity we found better agreement 
between the analytic and numerical results. 
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Figure 5: (Color online) Numerical spectra of coherently scat- 
tered light obtained by neglecting the third state in the model 
shown in Fig[l] (a) the spectrum from the two-level hydrogen 
atom and the strength-frequency parameters as in Fig[21 (b) 
the spectrum from the two-level Aj'*' ion and the strength- 
frequency parameters as in Fig[21 
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